Constrained one-step material decomposition reconstruction of head CT data from a silicon photon-counting prototype.

BACKGROUND: Spectral CT material decomposition provides quantitative information but is challenged by the instability of the inversion into basis materials. We have previously proposed the constrained One-Step Spectral CT Image Reconstruction (cOSSCIR) algorithm to stabilize the material decomposition inversion by directly estimating basis material images from spectral CT data. cOSSCIR was previously investigated on phantom data. PURPOSE: This study investigates the performance of cOSSCIR using head CT datasets acquired on a clinical photon-counting CT (PCCT) prototype. This is the first investigation of cOSSCIR for large-scale, anatomically complex, clinical PCCT data. The cOSSCIR decomposition is preceded by a spectrum estimation and nonlinear counts correction calibration step to address nonideal detector effects. METHODS: Head CT data were acquired on an early prototype clinical PCCT system using an edge-on silicon detector with eight energy bins. Calibration data of a step wedge phantom were also acquired and used to train a spectral model to account for the source spectrum and detector spectral response, and also to train a nonlinear counts correction model to account for pulse pileup effects. The cOSSCIR algorithm optimized the bone and adipose basis images directly from the photon counts data, while placing a grouped total variation (TV) constraint on the basis images. For comparison, basis images were also reconstructed by a two-step projection-domain approach of Maximum Likelihood Estimation (MLE) for decomposing basis sinograms, followed by filtered backprojection (MLE + FBP) or a TV minimization algorithm (MLE + TVmin ) to reconstruct basis images. We hypothesize that the cOSSCIR approach will provide a more stable inversion into basis images compared to two-step approaches. To investigate this hypothesis, the noise standard deviation in bone and soft-tissue regions of interest (ROIs) in the reconstructed images were compared between cOSSCIR and the two-step methods for a range of regularization constraint settings. RESULTS: cOSSCIR reduced the noise standard deviation in the basis images by a factor of two to six compared to that of MLE + TVmin , when both algorithms were constrained to produce images with the same TV. The cOSSCIR images demonstrated qualitatively improved spatial resolution and depiction of fine anatomical detail. The MLE + TVmin algorithm resulted in lower noise standard deviation than cOSSCIR for the virtual monoenergetic images (VMIs) at higher energy levels and constraint settings, while the cOSSCIR VMIs resulted in lower noise standard deviation at lower energy levels and overall higher qualitative spatial resolution. There were no statistically significant differences in the mean values within the bone region of images reconstructed by the studied algorithms. There were statistically significant differences in the mean values within the soft-tissue region of the reconstructed images, with cOSSCIR producing mean values closer to the expected values. CONCLUSIONS: The cOSSCIR algorithm, combined with our previously proposed spectral model estimation and nonlinear counts correction method, successfully estimated bone and adipose basis images from high resolution, large-scale patient data from a clinical PCCT prototype. The cOSSCIR basis images were able to depict fine anatomical details with a factor of two to six reduction in noise standard deviation compared to that of the MLE + TVmin two-step approach.

Photon-counting x-ray detectors for CT.

The introduction of photon-counting detectors is expected to be the next major breakthrough in clinical x-ray computed tomography (CT). During the last decade, there has been considerable research activity in the field of photon-counting CT, in terms of both hardware development and theoretical understanding of the factors affecting image quality. In this article, we review the recent progress in this field with the intent of highlighting the relationship between detector design considerations and the resulting image quality. We discuss detector design choices such as converter material, pixel size, and readout electronics design, and then elucidate their impact on detector performance in terms of dose efficiency, spatial resolution, and energy resolution. Furthermore, we give an overview of data processing, reconstruction methods and metrics of imaging performance; outline clinical applications; and discuss potential future developments.

Silicon photon-counting detector for full-field CT using an ASIC with adjustable shaping time.

Purpose: Photon-counting silicon strip detectors are attracting interest for use in next-generation CT scanners. For CT detectors in a clinical environment, it is desirable to have a low power consumption. However, decreasing the power consumption leads to higher noise. This is particularly detrimental for silicon detectors, which require a low noise floor to obtain a good dose efficiency. The increase in noise can be mitigated using a longer shaping time in the readout electronics. This also results in longer pulses, which requires an increased deadtime, thereby degrading the count-rate performance. However, as the photon flux varies greatly during a typical CT scan, not all projection lines require a high count-rate capability. We propose adjusting the shaping time to counteract the increased noise that results from decreasing the power consumption. Approach: To show the potential of increasing the shaping time to decrease the noise level, synchrotron measurements were performed using a detector prototype with two shaping time settings. From the measurements, a simulation model was developed and used to predict the performance of a future channel design. Results: Based on the synchrotron measurements, we show that increasing the shaping time from 28.1 to 39.4 ns decreases the noise and increases the signal-to-noise ratio with 6.5% at low count rates. With the developed simulation model, we predict that a 50% decrease in power can be attained in a proposed future detector design by increasing the shaping time with a factor of 1.875. Conclusion: Our results show that the shaping time can be an important tool to adapt the pulse length and noise level to the photon flux and thereby optimize the dose efficiency of photon-counting silicon detectors.

Feasibility of unconstrained three-material decomposition: imaging an excised human heart using a prototype silicon photon-counting CT detector.

RATIONALE AND OBJECTIVES: The purpose of this study was to evaluate the feasibility of unconstrained three-material decomposition in a human tissue specimen containing iodinated contrast agent, using an experimental multi-bin photon-counting silicon detector. It was further to evaluate potential added clinical value compared to a 1st-generation state-of-the-art dual-energy computed tomography system. MATERIALS AND METHODS: A prototype photon-counting silicon detector in a bench-top setup for x-ray tomographic imaging was calibrated using a multi-material calibration phantom. A heart with calcified plaque was obtained from a deceased patient, and the coronary arteries were injected with an iodinated contrast agent mixed with gelatin. The heart was imaged in the experimental setup and on a 1st-generation state-of-the-art dual-energy computed tomography system. Projection-based three-material decomposition without any constraints was performed with the photon-counting detector data, and the resulting images were compared with those obtained from the dual-energy system. RESULTS: The photon-counting detector images show better separation of iodine and calcium compared to the dual-energy images. Additional experiments confirmed that unbiased estimates of soft tissue, calcium, and iodine could be achieved without any constraints. CONCLUSION: The proposed experimental system could provide added clinical value compared to current dual-energy systems for imaging tasks where mix-up of iodine and calcium is an issue, and the anatomy is sufficiently small to allow iodine to be differentiated from calcium. Considering its previously shown count rate capability, these results show promise for future integration of this detector in a clinical CT scanner. KEY POINTS: • Spectral photon-counting detectors can solve some of the fundamental problems with conventional single-energy CT. • Dual-energy methods can be used to differentiate iodine and calcium, but to do so must rely on constraints, since solving for three unknowns with only two measurements is not possible. Photon-counting detectors can improve upon these methods by allowing unconstrained three-material decomposition. • A prototype photon-counting silicon detector with high count rate capability allows performing unconstrained three-material decomposition and qualitatively shows better differentiation of iodine and calcium than dual-energy CT.

Resolution characterization of a silicon-based, photon-counting computed tomography prototype capable of patient scanning.

Photon-counting detectors are expected to bring a range of improvements to patient imaging with x-ray computed tomography (CT). One is higher spatial resolution. We demonstrate the resolution obtained using a commercial CT scanner where the original energy-integrating detector has been replaced by a single-slice, silicon-based, photon-counting detector. This prototype constitutes the first full-field-of-view silicon-based CT scanner capable of patient scanning. First, the pixel response function and focal spot profile are measured and, combining the two, the system modulation transfer function is calculated. Second, the prototype is used to scan a resolution phantom and a skull phantom. The resolution images are compared to images from a state-of-the-art CT scanner. The comparison shows that for the prototype 19 lp / cm are detectable with the same clarity as 14 lp / cm on the reference scanner at equal dose and reconstruction grid, with more line pairs visible with increasing dose and decreasing image pixel size. The high spatial resolution remains evident in the anatomy of the skull phantom and is comparable to that of other photon-counting CT prototypes present in the literature. We conclude that the deep silicon-based detector used in our study could provide improved spatial resolution in patient imaging without increasing the x-ray dose.

Digital count summing vs analog charge summing for photon counting detectors: A performance simulation study.

PURPOSE: Charge sharing is a significant problem for CdTe-based photon counting detectors (PCDs) and can cause high-energy photons to be misclassified as one or more low-energy events. Charge sharing is especially problematic in PCDs for CT because the high flux necessitates small pixels, which increase the magnitude of charge sharing. Analog charge summing (ACS) is a powerful solution to reduce spectral distortion arising from charge sharing but may be difficult to implement. We investigate correction of the signal after digitization by the comparator ("digital count summing"), which is only able to correct a subset of charge sharing events but may have implementation advantages. We compare and quantify the relative performance of digital and ACS in simulations. METHODS: Transport of photons in CdTe was modeled using Monte Carlo simulations. Energy deposited in the CdTe substrate was converted to electrical charges of a predetermined shape, and all charges within the detector pixel are assumed to be perfectly collected. In ACS, the maximum charge received over any 2 × 2 block of pixels was grouped together prior to digitization. In digital count summing (DCS), the charge was digitized in each pixel, and subsequently, adjacent pixels that detected events grouped their charge to record a single, higher energy event. All simulations were performed at the limit of low flux (no pileup). The default tube voltage was 120 kVp, object thickness was 20 cm of water, pixel pitch was 250 µm, and charge cloud modeled as a Gaussian with σ = 40 µm. Variation of these parameters was examined in a sensitivity analysis. RESULTS: Detectors that used no correction, DCS, and ACS misclassified 51%, 39%, and 15% of incident photons, respectively. For iodine basis material imaging, DCS exhibited 100% greater dose efficiency compared to uncorrected, and ACS exhibited an additional 111% greater dose efficiency compared to digital charge summing. For a nonspectral task, the dose efficiency improvement as estimated by improvement of zero-frequency detective quantum efficiency, DQE(0) was 10% for DCS compared to uncorrected and 10% for ACS compared to DCS. A sensitivity analysis showed that DCS generally achieved half the benefit of ACS over a range of conditions, although the benefit was markedly less if the charge cloud was instead modeled as a small sphere. CONCLUSIONS: Summing of counts after digitization may be a simpler alternative to summing of charge prior to digitization due to the relative complexity of analog circuit design. Over most conditions studied, it provides roughly half the benefit of ACS and may offer certain implementation advantages.

Count statistics of nonparalyzable photon-counting detectors with nonzero pulse length.

PURPOSE: Photon-counting detectors are expected to be the next big step in the development of medical computed tomography (CT). Accurate modeling of the behavior of photon-counting detectors in both low and high count rate regimes is important for accurate image reconstruction and detector performance evaluations. The commonly used ideal nonparalyzable (delta pulse) model is built on crude assumptions that make it unsuitable for predicting the behavior of photon-counting detectors at high count rates. The aim of this work is to present an analytical count statistics model that better describes the behavior of photon-counting detectors with nonzero pulse length. METHODS: An analytical statistical count distribution model for nonparalyzable detectors with nonzero pulse length is derived using tools from statistical analysis. To validate the model, a nonparalyzable photon-counting detector is simulated using Monte Carlo methods and compared against. Image performance metrics are computed using the Fisher information metric and a comparison between the proposed model, approximations of the proposed model, and those made by the ideal nonparalyzable model is presented and analyzed. RESULTS: It is shown that the presented model agrees well with the results from the Monte Carlo simulation and is stable for varying x-ray beam qualities. It is also shown that a simple Gaussian approximation of the distribution can be used to accurately model the behavior and performance of nonparalyzable detectors with nonzero pulse length. Furthermore, the comparison of performance metrics show that the proposed model predicts a very different behavior than the ideal nonparalyzable detector model, suggesting that the proposed model can fill an important gap in the understanding of pileup effects. CONCLUSIONS: An analytical model for the count statistics of a nonparalyzable photon-counting detector with nonzero pulse length is presented. The model agrees well with results obtained from Monte Carlo simulations and can be used to improve, speed up and simplify modeling of photon-counting detectors.

Compression of CT sinogram data by decimation in the view direction.

PUPOSE: In clinical computed tomography (CT), the image data is acquired during continuous rotation. If the time during which the signal is integrated (the frame time) is too long, the data is blurred in the view direction (i.e., azimuthal blur). This can be overcome by having a high angular sampling rate, but for systems with limited bandwidth, the increased amount of data can be a problem. In this paper, we evaluate the benefit of maintaining a high angular sampling rate on the CT gantry and performing a decimation (digital low-pass filtration followed by a downsampling) in the view direction before the bottleneck of the data transfer chain. METHODS: A theoretical evaluation of the effects of the decimation is presented and the implementation of the digital filter is discussed. The compression scheme is evaluated on image data of a CATPHAN® 504 phantom and a human skull phantom. RESULTS: It is shown that digital decimation can be used to compress data before read-out with more remaining data fidelity compared to having longer frame times. Specifically, the method is shown to preserve the detail in the reconstruction of the CATPHAN resolution patterns and the human skull phantom. It is also demonstrated that the method can be used to prevent aliasing artifacts. CONCLUSIONS: Decimation in the view direction is presented as an alternative to increasing the frame time for CT systems with limited bandwidth of the data read-out. The method can be used to either remove aliasing artifacts or preserve spatial resolution. The proposed compression scheme can be implemented on the CT gantry and thus reduce the bandwidth requirements on the data transfer.

A method for geometric calibration of edge-on detectors in a CT-gantry.

PURPOSE: Photon-counting edge-on detectors are currently being considered for use in clinical computed tomography (CT) systems. A method for geometric calibration of edge-on detectors mounted on a CT-gantry has been developed and evaluated. The method is complementary to the geometrical calibration methods developed for CT systems using flat-panel detectors and takes the extra dimension of the edge-on detectors (along the direction of the x-rays) into account. METHODS: The method uses projection images of a simple phantom together with geometrical arguments to accurately estimate the orientation and relative position of the edge-on detectors. Both computer simulations and experimental measurements were used to verify the method. RESULTS: It is experimentally demonstrated that the method can determine the orientation of the detector with an accuracy of 0.08°. The method is also shown to be insensitive to errors in the modeled parameters used in the algorithm. CONCLUSIONS: The presented method can accurately determine the orientation and relative position of edge-on detectors mounted on a CT-gantry and can be used to evaluate the detector mounting and to produce an accurate forward model of the imaging system. Also, the method has potential to reduce the dimensionality of the geometric calibration of the full CT system since the direction of the x-rays with respect to the detector is measured.

Angular oversampling with temporally offset layers on multilayer detectors in computed tomography.

PURPOSE: Today's computed tomography (CT) scanners operate at an increasingly high rotation speed in order to reduce motion artifacts and to fulfill the requirements of dynamic acquisition, e.g., perfusion and cardiac imaging, with lower angular sampling rate as a consequence. In this paper, a simple method for obtaining angular oversampling when using multilayer detectors in continuous rotation CT is presented. METHODS: By introducing temporal offsets between the measurement periods of the different layers on a multilayer detector, the angular sampling rate can be increased by a factor equal to the number of layers on the detector. The increased angular sampling rate reduces the risk of producing aliasing artifacts in the image. A simulation of a detector with two layers is performed to prove the concept. RESULTS: The simulation study shows that aliasing artifacts from insufficient angular sampling are reduced by the proposed method. Specifically, when imaging a single point blurred by a 2D Gaussian kernel, the method is shown to reduce the strength of the aliasing artifacts by approximately an order of magnitude. CONCLUSIONS: The presented oversampling method is easy to implement in today's multilayer detectors and has the potential to reduce aliasing artifacts in the reconstructed images.

Improved signal-to-noise ratio for non-perpendicular detection angles in x-ray fluorescence computed tomography (XFCT).

The standard imaging setup in x-ray fluorescence computed tomography detects the fluorescence emission at a right angle with respect to the axis of the excitation beam. In this paper we have studied how the detection angle affects the signal-to-noise ratio (S/N), which is a major factor influencing the low-contrast sensitivity of the imaging system. This is done for an imaging setup using a collimated detector and a pencil beam of excitation x-rays. An ideal detection process is simulated for a generalized imaging case with gold/platinum tracers and experimental measurements are performed using a diagnostic x-ray tube. For monochromatic excitation, the results indicate that order-of-magnitude improvements of the S/N can be achieved by optimizing the detection angle. The maximal S/N, when exciting with an energy just above the K-edge, is achieved for large detection angles, i.e. with the detector close to the source. The improvements also transfer to polychromatic excitation sources and the experimental results show up to four-fold improvements of the S/N when changing the detection angle from 90° to 150°. Also, the changes of the S/N behavior when switching the fluorescent tracer is briefly demonstrated. These results suggest that the choice of detection angle should be taken seriously in the design of future XFCT imaging systems.